Multiloop calculations with Implicit Regularization in massiess theories 



E. W. Dias^'^E A. P. Baeta ScarpelliW L. C. T. Brito^l and H. G. Fargnoli( a H 

(a) Federal University of Minas Gerais - Physics Department - ICEx 
P.O. BOX 702, 30.161-970, Belo Horizonte MG - Brazil 
(b) Centro Federal de Educacao Tecnologica - MG 
Avenida Amazonas, 7675 - 30510-000 - Nova Gameleira - Belo Horizonte -MG - Brazil 

(c) Universidade Federal de Lavras - Departaraento de Ciencias Exatas 
Caixa Postal 37, 37.200-000, Lavras, Minas Gerais, Brazil and 

(d) Universidade Federal de Sao Joao del Rei - Campus Alto Paraopeba 
CAP Rod. MGU3, km 7 - 36420-000 - Ouro Branco/MG - Brazil 

(Dated: June 24, 2010) 

We establish a systematic way to calculate multiloop amplitudes of infrared safe massiess models 
with Implicit Regularization (IR), with a direct cancelation of the fictitious mass introduced by the 
procedure. The ultraviolet content of such amplitudes have a simple structure and its separation 
permits the identification of all the potential symmetry violating terms, the surface terms. Moreover, 
we develop a technique for the calculation of an important kind of finite multiloop integral which 
seems particularly convenient to use Feynman parametrization. Finally, we discuss the Implicit 
Regularization of infrared divergent amplitudes, showing with an example how it can be dealt with 
an analogous procedure in the coordinate space. 

PACS numbers: ll.10.Gh, 11.15.Bt, 11.25.Db 



I. INTRODUCTION 

Higher order calculations in massiess theories can be performed with the help of several techniques 
[H-[lj[j. In this paper we show that IR is a good tool for treating typical multiloop massiess integrals. 
Among the reasons that justify the use of IR, we can cite two important ones. The first one is the fact that 
the method works in the physical dimension of the theory, and this avoids complications with theories 



which are sensible to dimensional modifications (llj-[l3|. The second one is the simple algorithm that 



the IR provides for the identification of potentially symmetry violating terms: the surface terms which 
come from basic divergent integrals with Lorentz indices [lfl, [lH, [2(|, [24|-[26l|. 

It is just the second reason referred above that could bring some difficulties when one intends to apply 
the method to massiess models. This is because the IR standard expansion, used to separate the divergent 
from the finite part of an integral, in this case can only be performed with the introduction a fictitious 
mass; it turns out that the two parts are infrared divergent (here we are talking about infrared safe 
integrals). Of course, the parts must be added in order to cancel this mass; this is accomplished by 
means of a scale relation which introduces an independent mass parameter. However, the basic divergent 
integrals are a simple form of expressing the ultraviolet divergent content of the amplitude, since it does 
not ask for an explicit regularization [30(. 

On the other hand, this enforces the necessity of a method for performing the calculation of some 
pieces of the finite part obtained by such expansion, or the cancelation of the fictitious mass will not be 
evident. In this paper, we establish a systematic way to calculate multiloop integrals in massiess models. 
In addition, we develop a technique for dealing with an important kind of ultraviolet finite integral which 
emerge from the expansion of the integrand. We use a simple algebraic identity to put the integrand in 
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a convenient way to use Feynman parametrization. As a consequence, we obtain a direct cancelation of 
the fictitious mass when the finite and divergent parts are put together. 

It is well known that the use of a fictitious mass is not consistent when a genuine infrared divergence is 
to be treated. For the sake of completeness, we briefly comment on a procedure to implicitly regularize 
such kind of amplitudes, as it will be presented in [3l|. 

This paper is organized as follows: in section two we overview the basics of Implicit Regularization; the 
systematic calculation of multiloop massless integrals is presented in section three; we give an example of 
a five-loop amplitude in section four; some comments on the Implicit Regularization of genuine infrared 
divergent amplitudes are presented in section five; the concluding comments are in section six. Detailed 
calculations of a typical finite part are presented in Appendix A and the procedure for obtaining the scale 
relations is presented in the appendix B. 



II. OVERVIEW OF IMPLICIT REGULARIZATION 



Now, some important aspects of IR are needed; the detailed discussion can be found in the references 
[I3|-[29j. The condition is that the regularization, which is maintained implicit, does not modify the 
dimension neither the integrand. So, a simple cutoff is a good choice, since the basic divergences will not 
be calculated, and IR has a simple recipe to enforce symmetries. 

For treating massless theories, the following steps should be carried out: 

• the space-time and internal group algebra are performed and the amplitude is written as a combi- 
nation of basic integrals; 

• taking into account that the integral is infrared safe, a fictitious mass /i 2 is introduced. This must 
be done because the characteristic expansion of Implicit Regularization breaks the integral in parts 
that are infrared divergent; 

• the expansion of the integrand is carried out until the divergent piece does not have the external 
momentum in the denominator. The following identity is recursively used: 

1 _ 1 Pi - 2 Pi ■ k 

{ Pi - k) 2 - /i 2 (k 2 - y?) (k 2 - n 2 ) [{ Pi - k) 2 — [i 2 ] ' [ ' 

with pi a linear combination of the external momenta; 

• the divergent part is written in terms of basic divergent integrals of the type 

h °° (fI >-J k (fc 2 -,iY ln 1, a 2 y [ ' 

with 7 = 2/3 — 4, /, = J d 4 k / (2tt) 4: and where the A is to indicate the presence of an implicit 
regularization; 

• the basic divergent integrals with Lorentz indices are expressed as functions of basic divergent 
integrals without indices and surface terms; 

• the finite part is calculated; 

• scale relations are used in order to write the basic divergent integrals in terms of a non-null, arbitrary 
mass scale A 2 . The scale relation will allow an interplay between the finite and the divergent parts; 



• the limit /z 2 — » is taken. 
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III. THE CALCULATION OF MULTILOOP INTEGRALS 



A complete discussion on the ultraviolet content of infrared safe massless models within the context 
of Implicit Regularization is carried out in ref. (30j . It has been shown that the divergent part of an 
amplitude for such a theory can be completely displayed in terms of //*),(A 2 ), according to the definition 



(fc 2 - /i 2 ) ; 



■In* 



A 2 



(3) 



It has been also shown that after judiciously applying the identity ([T]), the divergent part of any multiloop 
amplitude will be contained in a set of integrals of the type 



Ju 



kll 



A 2 



(4) 



with p\ some linear combination of the external momenta and 2a < r + 3. In the expression above, it will 
be necessary the introduction of a fictitious mass /x 2 . A fictitious mass may always be introduced as long 
as the integral is infrared safe. This is necessary when using IR to treat massless infrared safe ultraviolet 
divergent integrals, because the expansion of the integrand, as we will see below, breaks it in two infrared 
divergent pieces. When a genuine infrared divergence appears this procedure can be problematic. For 
such cases a new procedure within IR defining basic infrared divergent integrals is necessary [3l| , and we 
will comment on this approach in section V. 

There are other finite contributions to the finite part besides the ones that come from the integral 
of eq. (H]). For some of them there is no analytic solution, but they have no problems with the limit 
/i 2 — > 0. This fact can be simply showed by proving that Q is infrared finite, since the whole amplitude is 
infrared finite. So, we are interested here in the explicit calculation only of this part of the finite content, 
because it is this one that will guaranty us, when the ultraviolet divergent part is considered together, the 
infrared finiteness of the amplitude. This justifies and shows the consistency of IR for treating infrared 
safe massless models. As discussed in the paper [30], the calculation of higher order renormalization 

group functions can be completely carried out by knowing the coefficients of the ^/og(A 2 )'s which display 
the ultraviolet content of all infrared safe massless amplitudes. Nevertheless, we must be secure that no 
new problem will emerge when the separation of these objects is performed. This is the reason why we 
dedicate ourselves to establish a procedure for solving the finite part of integrals of the type (@| . 
We now perform a complete calculation of a typical n + 1-loop integral, which has the form of (U): 




Carrying out the expansion of the integrand: 



ln n 



fc 2 — fi 2 
A 2 



Ik (fc 2 -M 2 ) 

=2A ( ;;i(M 2 )+^ +i) . 

The last term, the finite part, is given by 



H 2 )[(p-k) 2 



p 2 — 2p ■ k 



■ In" - 



k 2 - /i 2 



(5) 



(p 2 - 2p ■ fc) 2 



(fc 2 - M 2 ) (fc 2 -^ 2 ) 2 (fc 2 - /i 2 ) 2 [(p - k) 2 - ,i 2 



k a (p 2 - 2p ■ kf 



ln r 



fc 2 - fi 2 



l k (k 2 ~n 2 ) 3 l(p-k) 2 -v 2 } V A 2 
We first turn our attention to the tensorial divergent integral, which is given by 



(G) 
(7) 



1 log /iiAP > ~ 



k (fc 2 ~M 2 ) 3 



+20' - 1) 



k (fc 2 - m 2 ) 3 



ln J - 



A 2 

fc 2 - M 2 
A 2 



_ 1 J f A 1 1 : ,-lf k 2 -H 2 

-~4} 9 ^J k (fc 2 -M 2 ) 2 I A 2 



dk» 



{k 2 -ii 2 ) 2 



ln J - 



fc 2 - /i 2 
A 2 



(8) 
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The procedure is recursively repeated for l\ % og ~} v , l\ % og ~^ etc, until we obtain 

'2 = 9 -f"L ^few^V) + surface terms - ^ 

Recall that we still have to deal with the fictitious mass, which in the limit /i 2 — > will give infrared 
divergent pieces both in the ultraviolet divergent and finite parts. This problem is simply dealt with 
by the use of regularization independent scale relations (they can be easily obtained with the help of a 
cutoff), which read 

/S(A 2 ) = /S(A 2 )^E^m fc (^), do) 

k=l ' ^ ' 

with b = i/(47r) 2 , for arbitrary non-vanishing A. For infrared safe models a systematic cancelation of all 

powers of In between the ultraviolet divergent and finite parts finally crowns A a renormalization 

group scale. 

We use equations ^ and (fTUj) to write 

r(n+l)/,2s 9a/3 -A 1 n! J {k+1) 2 iV^l '^U 
fc=0 L i=l ' v ' ) 

fe=0 fe=0 i=l V 7 

After some algebra and the reorganization of the summations, we obtain 

C> 2 ) = ^ E (^S^) - Kl - 2™) ± In (£) } • (12) 



We see in the equation above that the second part, which is ultraviolet finite, diverges in the limit /j, 2 — >• 0. 
This part must be canceled by some contribution coining from the ultraviolet finite integral. We now 
turn ourselves to this integral. There is an interesting trick that allows us to use the traditional Feynman 
parametrization for solving this integral. We use the identity 

lna = lim - (a e - 1) (13) 

to write 

7 (n+1) = lim — V(-l)"- fc I k (14) 

e ™ £ "^ k\(n-k)\ a ' 1 1 

with 

rk_ 1 f fc a ( P 2 -2p-fc) 2 

" (-V) ke Jk (fc 2 -M 2 ) 3 " fce [(p-fc) 2 -At 2 ]' 

We leave the detailed calculation of this finite part for the appendix. We select the term of order n in e, 
since it is the only one which contributes. The result is 

= ^ E { [1 - 2i ^ 2 ] J ln * (£) - ("D"-* 1 2^+4 ln * ("£) } W 
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Remembering that is given by 

and using 



j(*0 



fe=0 



fe!(n-fc)! a 



£<-*: 



(n-fc)_ 



fc!(n-fe)! 



(17) 



(18) 



we see the perfect cancelation of the dependence on the fictitious mass [i. The final result is given by 

- {g (^jJ^*) - i ta< (-£)) 

+b [l - (l + (-1)" +1 ) 2- (n+2) ] } . (19) 



The procedure adopted above can be applied to any integral of the type of (g]) . The evaluation of multiloop 
amplitudes in massless theories becomes simple, since the divergences and the finite piece necessary for 
the cancelation of the fictitious mass are always originated in such kind of integral. In the next section 
we show an example in which the only necessary result is given in equation (TH 



IV. A MULTILOOP AMPLITUDE 



We now perform, as an example, the calculation of a multiloop amplitude. Specifically, we will treat 
the n-loop nested selfenergy of spinorial QED, as represented in the figure[T] With the subtraction of the 
subdivergences, we can write 

where the tilde is to designate the finite part of the (n — l)th order graph. From the Lorentz structure, it 
is easy to see that such an amplitude, whatever the order, may be displayed as T,( l '(p) = ftS^(p 2 ), with 
S^ % \p 2 ) a scalar function of p 2 . So, we obtain 

S (B) (P) = -27°/ fc A k 2 { p a _ k? ^ 1] (k% (21) 
and again the tilde is used to designate the finite part. Now, let us suppose that 

S^(k 2 ) = J2^ (— ■ (22) 



i=0 



So, we will get 

n 

E<">(p) = -2^ (23) 



according to the definition of the previous section. Consequently, from eq. (|19[) . it is found 

^«)(fc 2 )=EMn l (-g). (24) 



G 




FIG. 1: n-loop nested selfenergy of spinorial QED. 



From a simple calculation for the one-loop graph, we have 



(4tt) 2 



(25) 



and by induction, it is proved that the result (f2"3"]l is correct. As a specific result, for n = 5 it is simple to 
obtain 



= -&V I— 1 {5) ~ + —I {3) - —I {2) + 

w ' \i2 a 6 4 4 4 J 



(26) 



V. COMMENT ON IMPLICIT REGULARIZATION OF GENUINE INFRARED 

DIVERGENT AMPLITUDES 



The Implicit Regularization of infrared divergences is the subject of a work in progress [31| . Here 
we present the basics of the whole process. The essential steps which rendered Implicit Regularization 
adequate in the case of ultraviolet divergences have their counterpart for infrared ones. Moreover, a new 
scale appears, typically an infrared scale which is completely independent of the ultraviolet one. 

First, let us consider the following ultraviolet divergent massless integral and its result within Implicit 
Regularization, 



I = 



1 



(2tt) 4 k 2 {p-k) 2 
and then let us proceed with the following calculation 

u ■ d ' k 



I log (y)-b\n[-^ ) +26, 



P 



(2?r) 4 fc 4 (p- fc) 2 



(27) 



(28) 



By power counting U is infrared divergent and ultraviolet finite. In order to be able to use all the 
mathematics developed for ultraviolet divergent integrals we firstly note the following: 



1 



A d 4 ; 



1 



(2tt) 4 z 2 {z - uf 



(29) 



where z and u are configuration variables. 

Note the strinking similarity between the above z integral and equation (|2"T|) . We can thus write the 
result immediately: 



I(u 2 



1 



z 2 (z — u) 2 
h og {\- 2 ) -Mn(-u 2 A 2 ) + 2b. 



(30) 



But now Ii og {\ 2 ) is an infrared basic divergent integral and a scale relation has been used in order to 
get rid of a fictitious length in favor of the infrared scale I 2 = 1/A 2 . 
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Using this result and (j2l?)) in ([25)1 we have 



1 ' A 



A fA pA e i(p~k)x 



(4tt) 2 
i 



<?u / rf 4 x " ^ 2 e lku I{u 2 ) 

A „ — ipu 



(4tt) 2 
-&lnf-u 2 A 2 ) +26 



^ ^ (W* ) 



-^^o 5 (A- 2 )+61n^-^j+26j, (31) 

with A 2 = jl^A 2 , where 7 = 0, 5772... is the Euler-Mascheroni constant. This framework is similar to the 
Differential Renormalization of infrared divergences [l^l • 



VI. CONCLUSION 



This paper has been devoted to systematize the Implicit Regularization of massless models for any loop 
order. For the infrared safe theories, a systematic way to evaluate multiloop Feynman integrals in the 
context of Implicit Regularization has been presented. The ultraviolet content of probability amplitudes 
have a simple structure and we can easily identify all the potential symmetry violating terms, the surface 
terms. We have obtained a general expression for the scale relations, which establish the interplay between 
the divergent and hnitc parts of the amplitude, playing a fundamental role in the elimination of terms 
which are dependent of the fictitious mass. In addition, we have developed a technique for evaluating 
an important kind of finite Feynman integral which are typical of higher order calculations with Implicit 
Regularization in massless theories. We have extend the usual Feynman parametrization for integrals 
which are not written in terms of rational functions of the momenta. Through a simple example, we have 
exhibited the main elements of the procedure, showing how the finite part and the scale relations work 
together in order to restore the infrared safety of the amplitudes. Finally, we have discussed the Implicit 
Regularization of infrared divergent amplitudes, showing with an example how it can be dealt with an 
analogous procedure in the coordinate space. 



VII. APPENDIX A - THE CALCULATION OF THE FINITE PART 

We carry out the detailed evaluation of I k , which can be written, after Feynman parametrization and 
integration in k, as 

lk a =(^( A + B + 0> (32) 



where 



A = j2^kT) Pai ~^ )ke ~ 2 jf dx [p2{l ~ 2x)]2x{l ~ x)2 ^ ((-7)) ' (33) 



B = „ u~Zn i^ PaC-M 2 )* 6 " 1 / dx(l- X f- kf X ( ) CM) 

and 



c = - (1 _ fce )( 2 _ fce) ^(-^- 1 Si ^ a - *r Vd - 2x) ] . (35) 
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In the equations above, we have b = i/(A-n) 2 and 

H 2 = p 2 x(l — x) — jj, 2 . 



(36) 



Before considering the limit fj 2 — > 0, some important manipulations have to be done. First, we observe 
that 

±H 2 =p 2 (l-2x). (37) 

So, some integrations by parts are performed until the exponent of H 2 is ke. After this we can write 
H 2 — > p 2 x(l — x) without problem. We will have, then, 



/* = -b 



1 




1 - ke_ 


Pa | 



1 



EL 

ke{2-ke) \\ 2 



f 

J a 



dx 



(1-*) 
ke 



+ (1 - 2x) 



p x 



(38) 



In this point, it is convenient to evaluate each term separately, and so, we label the first term by £ and 
the second one by Q (/£ = £ + Q. So, we can write 

.2 \ kt 



£ = b 



1 1 



1 



2ke 1 - ke 1 - M 



Pa 



(39) 



For small e, we can perform a binomial expansion in each term and after some algebra we obtain 

11 1 -,000c /Ur\3 

2~ke 



■l-kel-% ~ 2fce ^^ (fce) ' ( 2 ) 

2 j=0 i=0 

oo 



i=0 L 



(ke) 



i-l 



(40) 



(fee 



Now we expand the term I w ) , for e — > 0. After this, we have 



i=0 L v 7 



(fee)'" 1 



The coefficient of (fee)™ will be given by 



+ 



2 - 



1 

2 

n-l 



n+1 



+ 



1 + fee In 



^2,2 / A 



-(fee) 2 In 



A 2 



+ 



(41) 



1, 2 /V 

ln 'a 2 " 



1 , „ / M 2 



In" 



A 2 



+ 2^-^ T7 ln"+ 1 ^ 



(n + 1)! 

Finally, the term of £ of order n in e can be written as 

£(") = 6 Pa ^[l-(2r"- 2 ]iln i (^)(fcer. 



i=0 



For C, we have 



C = -bp 



dxdx 



1 - ke 



1-x 
ke 



+ (1 - 2.t) 



o \ fee 
T 2 " 



which, in the rr order can be expressed as 



n+l 



i=0 



(42) 



(43) 



(44) 



(45) 
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VIII. APPENDIX B - THE SCALE RELATIONS 



In this appendix we discuss the main steps which are necessary in order to obtain the scale relations 
at n th order. As we have seen before, although the integral ([5]) by itself is infrared finite, when the 
separation by means of the relation (TTJ) is performed, we are left with two infrared divergent parts. The 
scale relations are important because they establish the connection between the finite and divergent parts 
in order to make the limit /i 2 — > possible. 

A typical basic logarithmic divergence at (n + l) th order in massless theories can be written, using a 
cutoff regulator, as 

and, going to the Euclidean space, it can be rewritten as 

d x i n ( X \ ..2 1 f 1 in ( x \ 



The logarithmic scale relations can be obtained from the following difference 



1 ln" +1 ^ ^ 



n + 1 



"+1 i / 2 \ 



i=0 



Finally, we write 



n+l 



At one loop (n = 0), for example, we have: 



i\ \ A 2 

i=0 



/,,,,i/^W,,,,^i-fcli.(£j]. ( r )( i) 



In two loops (n = 1), we obtain: 



Acknowledgements 



A. P. Baeta Scarpelli is grateful to CNPq for the financial support. The authors thank Marcos Sampaio, 
M. C. Nemes and B. Hiller for relevant discussions on this work. 



[1] A. V. Kotikov, Phys. Lett. B 375, 240 (1996). 



10 



K. G. Chetyrkin, A. L. Kataev and F. V. Tkachov, Nucl. Phys. B 174, 345 (1980). 
S. Laporta, Int. J. Mod. Phys. A 15, 5087 (2000). 

D. I. Kazakov, Theor. Math. Phys. 62, 84 (1985) [Teor. Mat. Fiz. 62, 127 (1984)]. 

A. Aste, Lett. Math. Phys. 77, 209 (2006). 

S. G. Gorishnii and A. P. Isaev, Theor. Math. Phys. 62, 232 (1985) [Teor. Mat. Fiz. 62, 345 (1985)]. 
N. I. Usyukina, Phys. Lett. B 267, 382 (1991). 

K. G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B 192, 159 (1981). 
I. Gonzalez and I. Schmidt, Nucl. Phys. B 769, 124 (2007) 

B. M. Kastening and H. Kleinert, Phys. Lett. A 269, 50 (2000). 

G. 't Hooft and M. J. G. Veltman, Nucl. Phys. B 44, 189 (1972). 
W. Siegel, Phys. Lett. B 94, 37 (1980). 

D. Stockinger, JHEP 0503, 076 (2005). 

O. A. Battistel, PhD thesis, Federal University of Minas Gerais (2000). 

C. R. Pontes, A. P. Baeta Scarpelli, Marcos Sampaio and M. C. Nemes J. Phys. G 34, 2215 (2007). 

E. W. Dias, A. P. Baeta Scarpelli, L. C. T. Brito, M. Sampaio, M. C. Nemes Eur. Phys. J. C 55, 667-681(2008) 

C. R. Pontes, A. P. Baeta Scarpelli, Marcos Sampaio, J. L. Acebal, M. C. Nemes, Eur. Phys. J.C 53, 121-131 

(2008) . 

O. A. Battistel, A. L. Mota, M. C. Nemes Mod. Phys. Lett. A 13, 1597 (1998). 

A. P. Baeta Scarpelli, M. Sampaio and M. C. Nemes, Phys. Rev. D 63, 046004 (2001). 

A. P. Baeta Scarpelli, M. Sampaio, B. Hiller and M. C. Nemes, Phys. Rev. D 64, 046013 (2001). 

M. Sampaio, A. P. Baeta Scarpelli, B. Hiller, A. Brizola, M. C. Nemes and S. Gobira, Phys. Rev. D65, 

125023 (2002). 

S. R. Gobira and M. C. Nemes, Int. J. Theor. Phys. 42, 2765 (2003). 

D. Carneiro, A. P. Baeta Scarpelli, M. Sampaio and M. C. Nemes, JHEP 0312, 044 (2003). 

M. Sampaio, A. P. Baeta Scarpelli, J. E. Ottoni, M. C. Nemes, Int. J. Theor. Phys. 45, 436 (2006). 

L. A. M. Souza, Marcos Sampaio, M. C. Nemes, Phys. Lett. B 632, 717 (2006). 

J. E. Ottoni, A. P. Baeta Scarpelli, Marcos Sampaio, M. C. Nemes, Phys. Lett. B 642, 253 (2006). 

E. W. Dias, B. Hiller, A. L. Mota, M. C. Nemes, Marcos Sampaio, A. A. Osipov, Mod. Phys. Lett. A 21, 
339 (2006). 

Brigitte Hiller, A. L. Mota, M. C. Nemes, A. A. Osipov, Marcos Sampaio, Nucl. Phys. A 769, 53 (2006). 
Orimar Battistel and G. Dallabona, Eur. Phys. J. C 45, 721 (2006). 

L. C. T. Brito, H .G. Fargnoli, A. P. Baeta Scarpelli, Marcos Sampaio, M. C. Nemes, Phys. Lett. B 673, 220 

(2009) . 

H. G. Fargnoli, A. P. Baeta Scarpelli, L. C. T. Brito, B. Hiller, Marcos Sampaio, M. C. Nemes, A. Osipov, 
Ultraviolet and Infrared Divergences in Implicit Regularization: a Consistent Approach, arXiv:1001.1543 
Javier Mas, M. Perez- Victoria, Cesar Seijas, JHEP 0203, 049 (2002). 



